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Abstract. 

The dynamics of the surface and inner atmosphere of the red supergiant star Betel- 
geuse are the subject of numerous high angular resolution and spectroscopic studies. 
Here, we present three-telescope interferometric data obtained at 11.15 fj.m wavelength 
with the Berkeley Infrared Spatial Interferometer (ISI), that probe the stellar surface 
continuum. We find striking variability in the size, effective temperature, and degree 
of asymmetry of the star over the years 2006-2009. These results may indicate an 
evolving shell of optically thick material close to the stellar photosphere. 



1. Introduction 



The red supergiant star Betelgeuse (a Orionis) is a key target for high angular resolution 
imaging studie s at optical and infrared w avelengths because of its large angular diam- 
eter (-0.047". Michelson & Peasel ll92lh and high apparent magnitude {Mk = -4.38). 
In recent years, a variety of interferometric campaigns have revealed up to three largely 
unresolved spots on the stellar siirface in various visi b le, near-infrared an d mid-infrared 
bands (e.g. I Wilson et all 1 19971 : iHaubois et all l2009l : iTatebeet al.ll2007l) . These spots 
range in intensity from a few to a few tens of percent of the total stell ar flux. The 
spots appear to be transient, with stability timescales of less than 8 weeks dWilson et al. 



19971 ). The properties of these spots have led to a consensus that they indicate large- 



scale convection cells, as predicted by ISchwarzschild (1975 ). 

Numerous puzzles remain regarding the effective temperatures of the spots, and 
their possible roles in mass loss from Betelgeuse and in shapi ng the asymmetric inner 
atmosphere of the star. Imaging studies in the ultraviolet (Gilliland & Dupreelll996l) 
suggest the presence of a plume of material at chromospheric temperatures. Wide-band 
radio images of Betelgeuse on spatial resolutions of a few stellar radii (/?«) also show 
asym metry in largely neutral gas at temperatures of a few thousand Kelvin dLim et al. 



Il998h . A plume at 6/?* was also observed by Kervella et al. (2009), possibly associated 
with the same gaseous component observed by ILim et all (119981) . Near- and mid-IR 
spectro-interferometric data for Betelgeuse have been modeled by a variety of authors 
using a dense shell of molecular gas at temperatures of ~2000 K located above a 3600 K 
photosphere ( Perrin et al."2004'), mainly composed of H2O (iTsujil I2OO0I : iPerrin et all 
20071 : lOhnaka et al.n200 9). This shell, known as a 'MOLsphere', is generally thought 



to lie between 1.3/?* and 1.5/?*, and has been shown to be patchy (Ohnaka et al. 2009, 
Evide nce for MOLspheres is also seen in other stars (e.g. T Lep, ILe Bouquin et all 
20090 . Betelgeuse is further known to possess an extended dust envelope, with two 



I 
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shells at angular radii of approxim ately l"and 2" ( Danchi et alJll994l) . A third, newly 
formed component was found by Bester et al. at 0.1". The mechanisms that 

drive the formation of the dust and shape the dynamic regions close to the surface of 
Betelgeuse are not well known. 



2. Observational methods and results 

The ISI is a three-telescope interferometer w ith heterodyne detect ion systems at each 
telescope (Wishnow et al., these proceedings: IWishnow et alfcolOh . This enables mea- 



surements over a narrow bandwidth of ~ 5. 4 GHz centered on a given CO2 laser local 
oscillator wavelength. The 11.15 fim wavelength chosen for these measure ments corre- 
sponds to a region that is free of prominent stellar molecular lines dWeiner et al.ll2003ljb . 

Visibilities and closure phases at spatial frequencies between 20-37 SFtQ were 
obtained over timespans of a few weeks in each year between 2006 and 2009 (see Table 
1 for an observing log). The van Cittert-Zernicke theorem relates the complex cross- 
correlation between the detected electric fields at different telescopes - the complex 
visibility - to the two-dimensional Fourier transform of the source brightness distri- 
bution. The ISI, like all optical/IR interferometers, however measures the Michelson 
fringe visibility, which is the power at a given spatial frequency normalized by the total 
source power. The closure phase is the sum of the fringe phases measured in each of 
the three baselines. While individual visibility phases var y due to atmo spheric turbu- 



lence, the closure phase is independent of these variations ( Jemiisonlll958 ). The spatial 
frequency sampled by a telescope baseline at a given time is proportional to b/A, where 
b is the baseline length orthogonal to the source direction and A is the wavelength. By 
obtaining measurements approximately every five minutes, the ISI samples visibilities 
and closure phases over a range of spatial frequencies and position angles as the Earth 
rotates. The telescope configuration used for the measurements reported here was an 
equilateral triangle with ~35 m baseline lengths. 



Table 1. Observing log. 



Year 


Dates (UTC) 


2006 


8, 9, 10 Nov; 7 Dec 


2007 


14, 15, 16 Nov 


2008 


22, 23, 24 Sep 


2009 


7 Oct; 4,5,7,9, 18,20, 22 Nov 



The visibility and closure phase data for each epoch are shown in Figure 1. The 
data were calibrated using nightly observations of Aldebaran (a Tauri), which was as- 
sumed to have a diameter of 20 mas and zero closure phase. Data that were affected 
by systematic errors and poor atmospheric conditions were discarded. The data in each 
epoch were fit using a simple image model of a centered uniform stellar disk, a point 
source at an arbitrary location and a large, spherically-symmetric dust shell that did not 
contribute to the visibility at the sampled spatial frequencies. The fits used a weighted 



' 1 Spatial Frequency Unit (SFU) =1x10^ cycles/radian. 
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Figure 1. Visibilities (left column) and closure phases (right column) from the 
2006, 2007, 2008 and 2009 epochs. Data in red, green and blue (diamonds, triangles 
and squares respectively) denote telescope baselines 1-2, 2-3 and 3-1, respectively. 
The closure phases are plotted against the position angle of the 2-3 baseline East of 
North. The solid lines indicate the best-fit image models. 



least-squares minimization technique, with the weights calculated as the inverse squares 
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of the errors. Each basehne and the closure phase data were weighted equally. For a 
baseline sampling the two-dimensional spatial frequency plane at a point (u, v), the 
model complex visibility was 



2nr Vm^ + 

where A is the fraction of the total flux contributed by the stellar surface, Ji denotes 
a Bessel function of order unity, r is the stellar angular radius, P is the fraction of the 
total flux contributed by the point, and x and y are the angular coordinates of the point 
with respect to the uniform disk center The fit results are presented in Table 2, with the 
errors in the last decimal places given in parantheses. The total flux density, Fv, of the 
star and dust was also estimated for each epoch assuming a flux density of 615 Jy for 
Aldebaran (iMonnier et al.lll998h . 



Table 2. Fit results. 



Year 


A 


r (mas) 


P 


X (mas) 


y (mas) 


(10^ Jy) 


2006 
2007 
2008 
2009 


0.51(2) 
0.55(2) 
0.47(2) 
0.540(5) 


24.5(8) 
25.7(9) 
24.6(8) 
25.9(1) 


0.05(1) 
0.02(1) 
0.04(2) 
0.013(1) 


-2.42(1) 
-3.25(2) 
-14.30(3) 
-24.40(1) 


-23.70(4) 
16.80(8) 
12.75(3) 
13.92(1) 


4.2(2) 
4.1(3) 
3.8(5) 
4.4(1) 



Figure 2 shows the best-fit model ima ges for each epoch. The 2006 data were 
previously analyzed by Tatebe e t al.l ( 20071) . and despite slightly different techniques 



used here, the results closely match. We find point sources of varying intensity in 
each epoch, as well as significant changes in the uniform disk radius of the star. Limb 
darkening is estimated t o redu ce the apparent stellar diameter by less than 1% at our 



wavelength dBester et alJI 19961) . and hence was not accounted for. 



2.1. The effective temperature of the stellar surface 

Our absolute measurements of the total flux of Betelgeuse, combined with the fitted 
stellar size and fraction of the total flux, enable an estimate of the effective temperature 
of the stellar surface. We used a Planck function, integrated over a 5.4 GHz bandwidth 
centered on 11.15/im, to estimate the temperature of the stellar disk. These results 
are presented in Figure 3. We did not account for extinction caused by the interstellar 
medium or by intervening layers of dust. However, the interstellar medium and the 
dust together are expected to attenuate the true stellar flux by only ~3. 1% iDanchi et all 



19941) . We also did not account for emission from the dust along the line of sight to the 



stellar face, again because of the low optical depth. 



3. Discussion 

3.1. The source of continuum opacity in the mid-IR 

In order to correctly interpret our results, it is crucial to elucidate the nature of the 
continuum opacity probed in our 11.15 /im observations. This is not an easy task. The 
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Figure 2. Depictions of the best-fit image models for Betelgeuse during each of 
the 2006, 2007, 2008 and 2009 epochs. Each figure includes the fit parameters: the 
fraction of the total flux from the star, the stellar radius in arcseconds, and the fraction 
of the total flux from the point. The point sources have been given the uniform disk 
sizes that they would have if they represented regions at a temperature of 7200 K. 
Our upper limit on the point source diameter is 20 mas. 



'surface' of a star at a given wavelength can be considered to be the location where the 
optical depth is unity. The stellar photosphere can be thought of as a surface where the 
opacity is equivalent to the Rosseland rnean opa city, which is a frequency average of 
the monochromatic opacity. IWeiner et all (12003 ah concluded that the mid-IR continuum 
opacity for solar abundance gas closely matches the Rosseland opacity under the stellar 
surface conditions present in red supergiant stars. This leads to the implicit conclusion 
that our 11.15 jum observations directly probe the stellar photosphere. This conclusion 
is further supported by the fact that ISI observations in a wavelength band inclusive of 
a strong H2O molecular line revealed a larger apparent stellar diameter for the AGB 
starMira (o Ceti) than that measured at 11.15 nm. This was indic at ive of a warm water 
shell that appeared transparent at 11.15/im ( Weiner et alj|2003bl) . iPerrin et al] ( 20071) 
further show that H2O molecular lines do not contribute significantly to the 11.15/im 
opacity in their model for the Betelgeuse MOLsphere. We can tentatively conclude that 
the contribution of a continuum of water lines from a shell above the stellar surface to 
the 11.15 jjm opacity in Betelgeuse is not significant. 

The contribution of the MOLsphere to the 1 1.15 nm opacity throu gh means other 
than water lines cannot however be neglected. IVerhoelst et all (l2006l) suggested that 
amorphous alumina (AI2O3) dust in the process of formation in the MOLsphere might 
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Figure 3. Estimated effective temperatures of the stellar surface of Betelgeuse, 
neglecting the contribution of any asymmetry, measured over the 2006, 2007, 2008 
and 2009 epochs. 



contribute to the opacity. Furthermore, 'Perri n et aP ( 20071 ) find, through modeling of 
spectro-interferometric data in the mid-IR obtained using the VLTI/MIDI instrument, 
that alumina as well as gaseous SiO appear to contribute to the opacity of the MOL- 
sphere. At 11.15 fim wavelength, this model suggests an optical depth of ~2, with the 

opacity dominated by alumina. 

An extension of the idea s of Reid & Mentenl ( r997h for the radio opacity of red 
supergiant and AGB stars by Tateb e & Townes shows that electron-hydrogen 

collisions in thermally ionized inner atmospheres could contribute significantly to the 
mid-IR opacity. This opacity is proportional to the temperature, naively implying that 
if a star became hotter, its apparent size would increase. Our present results provide 
evidence in support of such a phenomenon: the smallest size measurement corresponds 
to the lowest effective surface temperature (see Table 2 and Figure 3), and the larger 
size measurements correspond to higher effective surface temperatures. 



3.2. The origin of the spots 

A clear link has been made between spots o bserved at optical/near -IR wavelengths 
and the giant convecti on cell hypothesis of of ISchwarzschildl (119751) . as modeled by 
Chiavassa et al.l (120101) . However, it is harder to consider the spots we observe in the 
mid-IR as directly caused by high-temperature regions on the ste llar photosphere re - 
lated to convection cells. The spots seen at 700 nm wavelength by IWilson et all (Il997h . 
which contributed ~20% of the total iiux, were attributed to regions with temperature 
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excesses of at least 600 K from a mean photospheric temperature of 3600 K. The upper 
limit on the spot diameters was 15 mas. In contrast, Figure 2 shows that temperatures 
of at least twice the putative photospheric temperature of 3600 K are in some cases re- 
quired to produce spots of similar size that are consistent wit h our mid-IR observat ions. 
Such temperatures are inconsistent with the predictions of ISchwarzschildl (Il975[) for 
temperature fluctuations of ± 1000 K on red giant surfaces. 



4. Conclusions and future work 

Our present mid-IR continuum results for Betelgeuse show variability in the apparent 
stellar diameter, location and strength of asymmetries modeled as point sources, and in 
the effective stellar surface temperature, all on timescales of a year. Intriguingly, we 
find low effective surface temperatures, ~1000K below those usually associated with 
the photosphere, and an approximate scaling of the apparent radius with temperature. 
These results suggest that our measurements are dominated by the behaviour of cool, 
optically thick material above the stellar photosphere. A possible interpretation of our 
temperature and size results is that we have witnessed changes in this shell. This inter- 
pretation could be extended to explain the systemati c decrease in the apparent 11.15 yum 
size of Betelgeuse reported by iTownes et al.1 (l2009h over the interval 1 993-2008. The 



dominant 1 1.15 //m opacity so urce in the shell could, as suggested by IVerhoelst et al. 



and iPerrin et al] (l2007h . be alumina dust. If electron-hydrogen collisions, how 



ever, dominate the opacity, the apparent size changes could be attributed to temperature 
changes in largely neutral gas in thermal equilibrium at our measured effective surface 
temperatures. 

Clearly, further work is required to match our observations to radiative transfer 
models involving layers of material above the stellar photosphere, fully including the 
various possible opacity sources. More detailed modeling of the closure phase data is 
also necessary. Our modeled off-center point sources cannot be spots of similar size 
to those observed at optical/near-IR wavelengths, because the required surface temper- 
atures are unrealistic. The observed asymmetries in our data could possibly indicate 
large-scale asymmetry in the observed surface. 
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